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ABSTRACT 

A general solution of the Jeans equations for oblate scale-free logarithmic potentials 
is given. This provides all possible second velocity moments that can hold up a stellar 
population of flattened scale-free density against the gravity field. A two-parameter 
subset of second moments for the self-consistent density of Binney 's model is examined 
in detail. These solutions have the desirable property that the observable dispersions 
in the radial and proper motions can be given explicitly. 

In the spherical limit, the potential of these models reduces to that of the singular 
isothermal sphere. The Jeans solutions for scale-free densities of arbitrary fiattening 
that can correspond to physical three-integral distribution functions are identified. 

The problem of finding distribution functions associated with the Jeans solutions 
in fiattened scale-free logarithmic potentials is then investigated for Binney's model. 
An approximate solution of the collisionless Boltzmann equation is found which pro- 
vides a third (partial) integral of good accuracy for thin and near-thin tube orbits. 
It is a modification of the total angular momentum. This enables the construction of 
many simple three-integral distribution functions. The kinematic properties of these 
approximate DFs are shown to agree with a subset of the Jeans solutions — which 
are thereby confirmed as good approximations to physical solutions. The observable 
properties of the models are discussed. 

Key words: stellar dynamics - galaxies: kinematics and dynamics - galaxies: struc- 
ture 



1 INTRODUCTION 

Research in stellar dynamics applied to galaxies has gen- 
erally followed one of two fairly distinct methods. In the 
first, the statistical method, phase space is divided into a 
six-dimensional grid of elementary cells and the stellar sys- 
tem is described by a six-dimensional distribution function 
(DF) for which the collisionless Boltzmann equation gives 
the continuity condition. In applying this method, the full 
DF is often replaced by its low-order moments, which are 
governed by the Jeans equations (e.g.. Bacon 1985; Fill- 
more 1986; Evans & Lynden-Bell 1991; Magorrian & Binney 
1994). 

In the second, the orhit method, phase space is divided 
into a three-dimensional grid of three-dimensional slices, 
each of which contains one orbit. The stellar system is then 
described by a three-dimensional distribution function over 
all orbits (e.g., Schwarzschild 1979; Richstone 1980, 1984). 

In the first method, the motion of a star is described 
only by its instantaneous location, velocity and acceleration. 
This is just enough to follow the star through an elementary 
cell, but not enough to determine the orbit family to which 
the star belongs, or whether the orbit might be stochastic. 
In contrast, the second method requires the determination 



of each orbit over its full extent in phase space, and thus pro- 
vides a view of the structure of phase space for the potential 
under consideration. 

The present paper applies exclusively the first, statisti- 
cal, method. Within this method, we use approaches based 
on both the Jeans equations as well as the Boltzmann equa- 
tion. The two approaches complement each other. The 
advantage of the Jeans equations is that they are directly 
related to observables, such as the shape and tilt of the ve- 
locity ellipsoid. But, the drawback is that such solutions are 
worthless unless they correspond to a physical (i.e., every- 
where positive) DF satisfying the Boltzmann equation. 

Numerical schemes are already under development (e.g., 
Dehnen & Gerhard 1993; Binney 1994) that may enable gen- 
eral DFs to be built for large classes of galaxy models. In 
this paper, our attention is concentrated on just one class - 
the oblate scale-free logarithmic potential (Richstone 1980; 
Binney 1981; Toomre 1982). This potential possesses the 
astrophysically important property of a flat rotation curve. 
The most widely-used member of the class — often called 
Binney's model — has spheroidal equipotentials and a sim- 
ple self-consistent two-integral DF f{E, L],) (Toomre 1982; 
Evans 1993). We shall show here how to build simple ap- 
proximate three-integral DFs as well. This supplements the 



stone (1985a, b), who found numerical DFs for this potential 
by the orbit method. 

Section 2 describes the Jeans approach, and gives the 
complete solution of the Jeans equations for axisymmetric 
scale-free logarithmic potentials, as well as a useful two- 
parameter subset for the special case of Binney's model. 
Section 3 specialises these results to the limit of a spher- 
ical potential, and uses the Boltzmann approach to con- 
struct exact distribution functions. Section 4 describes the 
Boltzmann approach for flattened potentials, and presents a 
partial third integral which has fairly good accuracy in the 
part of phase space occupied by thin and near-thin tubes. 
The results from the two approaches are compared in Sec- 
tion 5. Finally, Section 6 describes how the observables can 
be obtained from the Jeans results. 



2 THE JEANS APPROACH 

We will work mainly in standard spherical coordinates 
{r,9,(l>), with 9 measured from the axis of symmetry and 
(f) the azimuthal angle. 

2.1 Jeans equations for spheroidal models 

Consider an axisymmetric galaxy model with a gravita- 
tional potential $ = $(r, ^) and with a phase-space DF 
/ = f{i',(^,Vr,ve,v,f,). In equilibrium, / must satisfy the 
coUisionless Boltzmann equation (Binney & Tremaine 1987 
[hereafter, BT], p. 279, eq. 4P-2): 
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Multiplication by Vr, vg, and v^, respectively, and subse- 
quent integration over velocity space, gives the Jeans equa- 
tions. These relate the mass density associated with / 



p = J J J f '^^^ '^^'^ '^^'t'' 
and the elements of the stress tensor. 



(2.2) 



fvtVj Avr Ave Av^, 



to the forces 9$/9r and T-^d^ld9. 1 n an axisymmetric 
model (vrVcf,) = {vgv^) = by the symmetries of the indi- 
vidual orbits, and we are left with two non-trivial relations: 
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+ '^[^VrVe) + {{vl)-{vl))coi9]=-'^"-^. 

These two relations between the four stresses p{vl), p(vrVe), 
p{vg), p{v^) must be satisfied at any point in an axisymmet- 
ric equilibrium model with potential $ and mass density p. 



i\ soiuiion oi inese equations ior given p anu vf corresponus 
to a physical equilibrium model only if it is associated with 
a DF f(r,9,Vr,ve,v^) > 0. 

2.2 Logarithmic potentials and densities 

Models with scale-free logarithmic potentials of the form 



<!> = lvlln[r'g{9)], 



(2.5) 



play a special role in modern astrophysics (e.g., Toomre 
1982). This is because the circular velocity in the equatorial 
plane is constant and equal to Vc. For some applications, it 
is useful to consider tracer populations of stars moving in an 
external gravity field rather than the self-consistent stellar 
density implied through Poisson's equation. So, we take the 
density to have the general form: 
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(2.6) 



Here, 7 is a constant while g{9) and h{9) are arbitrary func- 
tions. Symmetry with respect to the equatorial plane re- 
quires g{w — 9) = g{9) and h{-K — 9) = h{9). We could 
incorporate the factor g^ [9] in the definition of h{9), but it 
is convenient for what follows to write p in the form (2.6). 

When p is the self-consistent density, then 7 = 2 and 
h(9) and g(9) are related by 



h(9) = g^ (9) + \g(9)g' (9) cot 9 

+ \g(9)g"(9)-\g'(9)g'(9), 



(2.7) 



where we have taken 47rG = 1 and Vc = 1- When 7 7^ 2 the 
density (2.6) will overwhelm the density associated with the 
potential (2.5) either at large radii (7 < 2) or at small radii 
(7 > 2), so that results for this case should be used with 
caution. 

The most widely-used logarithmic potential is the 
model studied by Toomre in the 1970s and subsequently 
popularised by Richstone (1980) and Binney (1981; see also 
BT, chap. 2). This has 



g{9) = sin^ 9 + 



cos' 9 



(2.8) 



so that the equipotentials are similar concentric spheroids 
with axis ratio q. The associated self-consistent density is 
of the form (2.6), with (BT, p. 46, eq. 2-54b) 



(2.3) h{9) = Q[2-g{9)] = Q[l + {1-Q) cos" 9], 



(2.9) 



where g(9) is given in (2.8), and we have written Q = ■ 
The density becomes negative when Q > 2. We shall refer 
to this potential and density as Binney's model. 

Other cases of interest with potentials of the form 
(2.5) include the singular isothermal spheroids, which have 
p = r~"(sin" 9 + Q cos' 9) and g{9) is the integral given in 
equation (B2) of Qian et al. (1995). A family of scale-free 
models with elementary g(9) and h(9) and nearly spheroidal 
density distribution is given in Appendix C of de Zeeuw & 
Pfenniger (1988). 

2.3 General solution of the Jeans equations 

The potential (2.5) and density (2.6) have the desirable at- 
tribute of scale-freeness. Their properties at radius r' = kr 
follow from those at radius r by a simple magnification — 



that are also scale-free, i.e., that satisfy f{kr, 9, Vr, vg, v^) = 
k"* f(r,0,vr,ve,v^) (e.g., Richstone 1980; White 1985). 
The associated stresses therefore have the following form: 
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where F\, F2, F3, and Fi are functions of 9. They must 
fulfill the following conditions 



Fi{7r-0) = Fi{0), 
F2(7r-0) = -F2(e), 



FsiTT-e) = Fsie), 

Fi{w-e) = Fi{e). 



(2.11) 



These constraints guarantee that the stresses are symmet- 
ric with respect to the equatorial plane. That is, p(v^), 
p(vg) and p(v^) must each have the same value at the point 
(R, z) = (r, 9) and its mirror image (R, —z) = (r, Tr—9). But, 
p{vrVg) must have opposite sign at two such mirror points 
because it defines the tilt of the velocity ellipsoid with re- 
spect to the radial direction in a meridional plane (R, z). 

Substitution of the forms (2.10) into the Jeans equa- 
tions (2.4) reduces them to two coupled first order ordinary 
differential equations: 

Fi{0)+ \cot0-%^]F2{e) 
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(2.12a) 



(2.12b) 



These equations provide two relations between the four func- 
tions Fi{9), . . . , Fi{9). We are free to pick two of the func- 
tions arbitrarily and solve (2.12) for the other two. We 
choose to prescribe Fi{9) and F2{0), subject to the con- 
straints (2.11) and the requirement Fi{9) > 0. This choice 
is motivated by the fact that it is then easy to consider the 
special case _F2 = 0, in which the velocity ellipsoid is every- 
where aligned exactly along the spherical coordinate system. 

We use equation (2.12a) to eliminate Fi{9) from 
(2.12b). It is straightforward to integrate the resulting first 
order differential equation for F3 . Substitution of the result 
in (2.12a) then also gives F^. We find: 



Fsie) = i{0) + J{0) + F2{0) cote, 

Fi{0) = h(0) - 1(0) - J(0) 

-{'i-2)F,{0)+g\0) 

with 



(2.13a) 
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tures for I{0) and J{0). It follows that 
Fs(0) = Fi(0) = iA(0) - 1(7 - 2)i^i(0), 



(2.16) 



so that (vg) = (v^) on the minor axis. Of course, this is re- 
quired by axisymmetry — on the minor axis the i/i-direction 
is identical to the ^-direction. This completes the derivation 
of the general solution of the Jeans equations for models 
with scale-free logarithmic potentials (2.5) and power-law 
densities of the form (2.6). It is particularly fortunate that 
the general solution can be written in simple analytical form. 

Inspection of equations (2.13) shows that I{0) and J{0) 
occur only in the combination I{0) + J{0), so that we could 
have combined them. However, the separate definitions of 
I{0) and J{0) shows the basic mathematical difference be- 
tween solutions with 7 = 2 and those with 7 7^ 2. In the 
former case J{0) vanishes, so that no integration of the ar- 
bitrary functions Fi{0) and F2{0) is required. As we will 
mainly be concerned with this case, we give the simpler form 
of the general 7 = 2 Jeans solutions explicitly: 



F3{0) = I{0) + F2{0)cot0, 

F^{0) = h{0)-I{0) + g'{0)- 



d0 [gH0) 



(2.17) 



where I{0) is again given in equation (2.14). In this case we 
have ^3(0) = Fi{0) = 1(0) = |/j(0). 



2.4 Orientation of the velocity ellipsoid 

The velocity ellipsoids defined by the general solution (2.13) 
are all aligned with the i/i-direction. But, in a plane of con- 
stant (f), they are misaligned from the r- and ^-direction 
when (vrVg) ^ 0, i.e., when F2{0) ^ 0. We define velocity 
components {v\,Vtj) in the meridional plane by 



cos sin Q \ Ivr 
- sin cos Q I \vg 



(2.18) 



where we choose at each position in such a way that 
(vxViy) = 0. This is accomplished by taking 

2_F2 



tan 20 
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Fi — F3 



(2.19) 



so that is a function of only, as expected in a scale-free 
model. The transformation (2.18) diagonalises the stress 
tensor. The principal components p(v\) and p(v^) are given 
by 



P{-"1) 



with 
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F±{0) = ^{Fs + Fi) ± ly'{Fs-Fiy +4FI 



(2.20) 



(2.21) 



A necessary — but not sufficient — condition for the stresses 
(2.10) to correspond to a physical equilibrium model is that 
FijFsjFi > 0, since they give the velocity average of the 
non-negative quantities v^, Vg, and v^. We must also require 
F- > and F^ > 0, which means we must have 



i^2 <FlF3, 

for the entire range of between and 7r/2. 



(2.22) 



We consider briefly the special case corresponding to a two- 
integral DF / = f{E,L],). In this model, the stellar veloci- 
ties have no preferred direction in the meridional plane. In 
other words, (yl) = (vg) and (vrVe) = 0, so that the velocity 
ellipse in the (R, 2)-plane is a circle. The associated stresses 
follow from the general solution by taking 



Fi{0) = Fs{0), 



F2{0) = 0. 



(2.23) 



Substitution of (2.23) in equations (2.13a) and (2.15) results 
in an integral equation for Fi which is of Volterra type. It 
is easily solved to give 



sin'' 61 
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(2.24) 



Substitution of this expression in equation (2.13b) then gives 
Fi(0) as: 



F,{0) = h{0) + {l-y)Fi{0). 



(2.25) 



When 7 = 2, Fi{0) = F3{0) = I{0), and Fi{0) = h{0)-I{0). 
These results can also be obtained from Hunter's (1977) so- 
lution of the Jeans equations for general f(E, ij) models. 

The DF of the scale-free two-integral models is of the 
form f(E,Ll) oc exp(-7_£:)/-,(»?^), where = L,/Lc(E), 
and Lc{E) is the value of Lz on the circular orbit of en- 
ergy E. The function can be evaluated by means of the 
contour method of Hunter & Qian (1993). A number of 
specific cases can be found in the literature. Toomre (1982) 
and Evans (1993) give the elementary result for the self- 
consistent Binney model. Qian et al. (1995) give f{E, i^) 
for scale-free spheroidal densities in Binney's potential, and 
also for the self-consistent scale-free spheroids. They show 
that the second moments for the latter models are elemen- 
tary when 7 = 2. 

2.6 A two-parameter solution for Binney's model 

For the specific case of Binney's model, the entire general so- 
lution of the Jeans equations can be written down explicitly. 
We restrict ourselves to oblate models with 1^2 < g < 1 
(i.e., 1 < Q < 2.). The self-consistent density has 7 = 2 so 
that J{0) = 0. Taking g{0) and h{0) as given in (2.8) and 
(2.9) respectively, the integral (2.14) is elementary. We find 
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h{0)-I{0) 



^0 + Q{2-Q) cos" 0], 
|[(2Q-1) sin" + Q{2-Q) cos" 0], 



(2.26) 



Substitution in equation (2.17) then gives F3{0) and Fi{0), 
for each choice of F2{0). Note that, when 7 = 2, Fi{0) does 
not even appear in the scale-free Jeans equations (2.12) and 
so its choice does not affect the other stresses. 

Picking (2.23) for Fi(0) and F2(0), and taking 7 = 2, 
we recover the stresses for the two-integral f{E, i") model 
derived by Evans (1993). They have a very simple functional 
dependence when written in terms of cylindrical coordinates 
{R, 4>, z), with R^ + = r" and z = r cos 0, namely 



, 9 , a,R + b, Rz + c, z 
p(v^) = i 



(2.27) 
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where , 6j and Cj are constants which depend on q ( j = 
R,4>,z). Evans & de Zeeuw (1994) showed that this simple 
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form (2.27) has a desirable property — it allows an explicit 
evaluation of the line-of-sight projected second moments, 
such as the dispersions in the radial velocities and the proper 
motions. We elect to study the subset of the general solu- 
tion (2.17) and (2.26) of the Jeans equations that possesses 
this pleasing characteristic. We now show that this defines 
a two-parameter family of anisotropic stresses for Binney's 
potential, and so this property is not restricted to the unique 
f{E, Ll) model. 

Appendix A gives the transformation between the 
stresses in spherical and cylindrical coordinates. It is 
straightforward to show that all Jeans solutions (2.17) and 



Fi{0) = 1(0) + H1+H2 g(0) sin" 0, 
F2{0) = H2g{0}sm0cos0, 



(2.28) 



lead to stresses of the desired form (2.27). Here, I{0) is 
given in equation (2.26) while Hi and H2 are constants. We 
emphasise that this is a very restricted subset of the general 
solution (2.17) — instead of two free functions Fi{0) and 
F2{0), we have merely two free parameters Hi and H2. So- 
lutions with Hi = correspond to velocity ellipsoids that 
are everywhere aligned with the cylindrical coordinate sys- 
tem. Solutions with H2 = are everywhere aligned with the 
spherical coordinate system. 

Substitution of the choice (2.28) in the general solution 
(2.17) yields 



F3{0) = I{0) + H2g{0)cos^0, 

Fi{0) = h(0) - 1(0) + _ff2(<?cos" 6i-sin" 0), 



(2.29) 



with I{0) and h{0) — 1(0) given in equation (2.26). It is help- 
ful to write this Jeans solution out explicitly in cylindrical 
coordinates to emphasise its simplicity: 
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where the relationships between the coefficients Aj , Bj , Cj 
and our two anisotropy parameters Hi and H2 are given in 
Table 1. 



riiiany, lu prepare rur ine curiiparrsun ur ine present 
results with those from the Boltzmann approach (Section 
5.1), we assemble our solutions for Binney's model for the 
case H2 = in the following form 
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These formulae will be useful later. 



2.7 Limits on the two-parameter solution 

The principal stresses p{v\) and p{v^) of the two-parameter 
Binney models described in the previous subsection are 
given in equation (2.20), where now 

F±{0) = I{0) + ^H^ + ^H2g{0) 



± y^Hf - 2H1H2 g(e) cos 261 + Hlg^(ff). 



(2.32) 



The tilt angle of the velocity ellipsoid in the meridional 
plane is given by 



tan 20 



H2 gjO) sin 20 
Hi-H2g{0)cos2e' 



(2.33) 



Necessary (but not sufficient) requirements for positive 
stresses are that p(v%}, p(v^) and p(v^) are non-negative. 
This requires that > 0, Ci > 0, C3 > and yl4 > 0, and 
translates into the following conditions on Hi and H2: 



-|Q(2-Q) < Hi, 

(|Q - 1) < < Q - I, 

-i < Hi + H2. 



(2.34) 



These define a region in the (_ffi , _ff2)-plane illustrated in 
Fig. 1. 

Necessary and sufficient conditions for positive stresses 
are that the principal components p(v\) and p(v^) are non- 
negative. In fact, it is slightly easier to work with the sum 
and product of the principal components. The sum equals 
p(v^) + p(v^), and hence is already positive, so only the 
product remains to be checked. This implies that at every 
point (r, 6,1/1), we need 



A(6) = a sin* + b sm + c>O, 
where 

a = {Ai - C\){As - Cs) + Bl 
h = C\{As - Cs) + Cs{Ai - C\) 

c = C1C3. 



B 



(2.35) 



(2.36) 



The requirement (2.35) already holds at 6 = and = ^ir, 
since -D(O) = C1C3 and D(^Tr) = A1A3. The Jeans solutions 
with negative stresses correspond to the cases when A(6) 
has a minimum in the range < 6 < ^-tt which is negative. 
Therefore, the additional regions that must be expunged 
from Fig. 1 are those simultaneously satisfying 



2H2=2q-2-l 




8H2=q-2-2 



Figure 1. The (/fi , /f2)-plane of two-parameter anisotropic solu- 
tions of the Jeans equations for Binney's oblate scale-free logarith- 
mic model with < g < 1. The solutions with non-negative 
stresses p(u^), p[v^) , and p{v^) lie in the hatched area. The re- 
quirement of positive principal stresses p{v^) and p{v^) leads to 
the exclusion of two tiny regions near the two corners in the third 
quadrant. These are invisible on the scale of the Figure. The 
filled circle indicates the stresses with Hi = H2 = 0, which are 
associated with the f[E,Lz) model. Spherically aligned velocity 
ellipsoids have H2 = 0, while cylindrically aligned velocity ellp- 
soids have Hi = 0. Approximate DFs for solutions with H2 = 
and Hi < are derived in Section 4. 



Each of these four inequalities defines an area in the 
(-ffi , -ff2)-plane. The boundary curves follow by substitut- 
ing the expressions in Table 1 in the definitions (2.36), and 
taking the equality signs in (2.37). This gives 
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H2 
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2(i-Q)2-4Hi' 
(l-Qf[Hi+Q(2-Q)] 
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H2 = {Q-l)Hi 
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a > 0, 



< -6 < 2a, 



-r +4:ac < 0. 



(2.37) 



The area that lies inside all four boundaries needs to be 
excluded from the region defined by (2.34). Somewhat sur- 
prisingly, this results in the deletion of only two tiny areas 
near the two corners of the hatched region that lie in the 
third quadrant of Fig. 1. We will ignore them because they 
lie away from the Jeans solutions with H2 = 0, for which we 
will construct approximate DFs in Section 4. 

The position of the two-integral solution (2.23) is 
marked in Fig. 1 — it lies at the origin in the {Hi,H2)- 
plane and corresponds to a positive definite DF (provided 
the flattening is not greater than q = ^\/2). As Hi 00, 



tended. It is clear from Fig. 1 that the set of Jeans solutions 
for Binney's self-consistent model is large. What is far from 
clear is which solutions are physically relevant. 




H.=0.5 



3 SPHERICAL POTENTIALS 

The existence of chaotic orbits prevents the extension of an 
analytic third integral to all the orbits in generic axisymmet- 
ric potentials. So, ascertaining which of the Jeans solutions 
have physical three-integral DFs is difficult in non-spherical 
models. We will consider this problem in Section 4, and 
here investigate the limiting case where the potential is ex- 
actly spherical. This allows the construction of galaxy mod- 
els with three-integral DFs and triaxial kinematics. This is 
possible because the total angular momentum L provides an 
exact third integral of motion in a spherical potential. 




3.1 General Jeans solutions 

In spherical symmetry, g(d) = 1, so all the potentials (2.5) 
reduce to that of the singular isothermal sphere (BT, p. 
228). The general solution for the stresses associated with 
a flattened density p = r~~' h{9) is of the form (2.10) with 
g{9) = 1. The functions Fi{9) and F2{0) are arbitrary, while 
F3(0) and Fi(0) follow from (2.13) upon taking g(0) = 1. 
However, it follows from the properties of the individual 
orbits in a spherical potential that {vrVg) = for each of 
them. So, all solutions with F2{0) ^ are unphysical! The 
remaining solutions have 



Fsie) 

F,{0) 
and 

1(0) = 
J(0) = 



I{0) + J{0), 

h{0)-I{0)-J{0)-{'(-2)F^{0), 



h(0)sin0cos0A0, 



(3.1) 



(2-7) 



(3.2) 







Fi{0)sm0cos0A0, 



where Fi{0), which defines the angular variation of the radial 
velocity dispersion, is arbitrary. Once it and h{0) have been 
chosen, the other stresses are fixed. In the special case 7 = 3 
we have (v^) = {vl) + {vl) + (vD = 1. 

Both h{0) and Fi{0) are even functions symmetric 
about = Tr/2 and so each may be expanded as a Fourier se- 
ries of the form ^ cos 2£^ (see, e.g., Mathews & Walker 
1964). These series can be rearranged as series in even pow- 
ers of sin 0, so we may write without loss of generality: 



h{0) 



Fi(0) 



(3.3) 



Figure 2. The (i^i , )-plane of two-parameter anisotropic so- 
lutions of the Jeans equations for the singular isothermal sphere, 
which is the limit g ^ 1 of Binney's model. The solutions with 
non-negative stresses lie in the hatched area. The filled circle in- 
dicates the stresses with Hi = H2 = 0, which are associated with 
the f{E) isotropic model. The only physical solutions are those 
with spherically aligned velocity ellipsoids, which lie along the 
thick solid line H2 = 0. The model built exclusively with circular 
orbits has H2 = and Hi = —1/2, while the radial orbit model 
has H2 = and Hi 00. 



White (1985) showed that the density r 
be reproduced by non-negative DFs of the form 



can 



7- 2m 7- 2n 

Vm,nF Fz exp 



where m + n > —1 and 2n > 
constant given by 



('/ + 2n + 2m)E], (3.4) 
1, and •r]m,n is a normalization 



(m + « + ^)'"+"+tr(« + l) 



7rT{n+^)T{n + m + l) 
The corresponding stresses are 
sin"" 



(3.5) 



1 



P{Vr 



2m + 2n + ~f ri 
m + n + 1 2x 

-P{vr), 



(3.6) 



ra -I- 1 
{2n + l)p{vl). 



For fixed 7 and ra, the stresses and the density have the 
same angular dependence for all m, but the anisotropy ra- 
tios {vg) I {vl) and (f|)/(fr) do depend on m. White showed 
that further DFs can be constructed from (3.4) by linear su- 
perposition, and that the general solution of the Boltzmann 
equation is 



f(E, F^,FI) = ^ am,nfrr 



(3.7) 



We now establish that it is possible to find DFs generat- 
ing the Jeans solution (3.1) that corresponds to each pair 
h{0), Fi{0). The specific case of scale-free spheroidal densi- 
ties is discussed in some detail by de Bruijne, van der Marel 
& de Zeeuw (1996). Kochanek (1994) has also recently stud- 
ied Jeans solutions of flattened denisties in the isothermal 
sphere. 



where the am,n are constants specifying the fraction con- 
tributed by each component. To recover the angular vari- 
ation h{0) of the density and Fi{0) of the radial velocity 
dispersion, as given in equation (3.3), we therefore require 



m = 



^ 2m-|-2ra-|-7 

m = 



bn, 



(3.8) 



iui n = u, 1, . . mere are iiiaiiy iiiuie uiiKiiuwiis cim,n 
than constraints a„,b„. So, there are many ways to define 
three-integral DFs that reproduce the required velocity dis- 
persions. Of course, we have not shown that these are all 
positive definite, although it is likely that at least some are - 
depending on the choice of h{9) and Fi{9). For example, de 
Bruijne et al. (1996) provide two families of three-integral 
DFs consistent with scale-free spheroidal densities of arbi- 
trary flattening, which cover a large range in anisotropy, and 
are all physical. 



3.2 A one-parameter solution for spherical models 

It is instructive to consider in some detail the limiting be- 
havior of the two-parameter Jeans solutions for the self- 
consistent Hinney model discussed in Section 2.6. The den- 
sity is now spherical, so that h(6) = 1, and 1(6) = i. Since 
7 = 2 we have J{0) = 0. We then find from equations (2.28) 
and (2.29): 



F2{e) = _ff2sin6icos6i, 
F3{0) = \ + H2Cos'e, 
Fi{0) = 1 _ff2(cos^ - sin^ 61). 



(3.9) 



These Jeans solutions provide anisotropic stresses for the 
singular isothermal sphere. They are positive in the region 
in the (_ffi , _ff2)-plane defined by (cf. eq. [2.34]) 



_i <Hi+H2 



(3.10) 



which is illustrated in Fig. 2, and may be compared with 
Fig. 1. However, as we have seen, all physical solutions must 
have F2 = because of the properties of individual orbits, 
so they all have H2 = 0. We are left with a one-parameter 
solution 



F,{0) = l + Hi, 
Fs{0) = h 



F2(0) = 0, 



(3.11) 



which is indicated as the thick solid line in Fig. 2. It follows 
from the derivation in Section 3.1 that each of these Jeans 
solutions corresponds to a DF of the general form (3.7) with 
ra = 0, so that (see also White 1985; Gerhard 1991; Kulessa 
& Lynden-Bell 1992) 

CO 

f(E, i") = am.oVm.oL^"" exp[-(2 + 2m)E], (3.12) 

m = 

with rjmfi given in equation (3.5), and we must have 

CO CO 

y«™.o = l, y ^ = l + 2H^. (3.13) 

^ ^ m-l-1 ^ ' 

m = Q m = Q 

Again, many DFs are possible. The simplest is obtained by 
taking only one component in the series (3.12): 

f{E, i') = »?„,oi"" exp[-{2 + 2m)E], (3.14) 

with 

2 Hi 



1 nis snows inai _oi is a measure oi me anisoiropy oi me 
model. Hi = — i corresponds to the circular orbit model 
(no radial velocity dispersion) and Hi = to the isotropic 
model, which has / = f{E). When Hi 00, the model 
approaches the radial orbit model. The DF (3.14) is non- 
negative and hence physical for each value of Hi in the al- 
lowed range. 



4 THE BOLTZMANN APPROACH 

We now investigate which of the two-parameter Jeans solu- 
tions can correspond to three-integral DFs for the flattened 
Binney models. 



4.1 Classical integrals 

The coUisionless Boltzmann equation (2.1) for the phase- 
space density / can be put into the form 



= VrA + veB, 
with 



(4.1) 



2 9$, 9/ 2 9/ 



df 



2^— + y''e + V4,-T — )— -Vg— -V^j;^, (4.2 



d 



dr dv 



and 

The derivatives of the potential (2.5) are (with Vc = 1) 

9$ _ 1 9$ _ g'{0) 

dr ~ r' 30 ~ 2g{0)' 

This axisymmetric potential has two classical integrals, the 
energy i? = $ -|- \{vr -\- Vg -\- v^) and the ^-component of 
the angular momentum = rv^sm0. As is often done for 
convenience sake, we use these two integrals in the following 
form 

— 2 E 1 



(4.3) 



(4.4) 



h 



r^g{0) 
I2 = L\ = v^tf, sin^ . 



(4.5) 



These two integrals each fulfill the Boltzmann equation (4.1) 
by giving A = Q and 5 = separately. 

4.2 A partial integral for Binney's model 

We now restrict ourselves to the specific case of Binney's 
model. Its potential does not have an exact third integral, 
as is evident from the existence of some stochastic orbits 
(e.g., Richstone 1982). Here, no attempt is made to find 
an approximate third integral of reasonable accuracy for all 
orbits. Instead, our aim is to construct an approximate in- 
tegral that has high accuracy for thin tube orbits. With 
such a partial integral, we can build models with tangen- 
tially anisotropic velocity distributions. 

In the spherical limit (g = 1), the potential has an 
exact third integral, namely the total angular momentum 
= r^ (vg + v^). A number of authors (e.g., Saaf 1968; In- 
nanen & Papp 1977) have suggested generalisations of i as 
approximate third integrals for nearly spherical potentials. 
Here, we introduce a new modification, namely 



I + 2H1' 



(3.15) J3=i2^(^^^^2(^2^^2^^(^^^ 



(4.6) 



laDie z. Accuracy oi tne partial integral is . ivleriaional cross sections ol 
the orbits are displayed in Fig. 3. 
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E 


0.9 .4800 .5000 
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.2969 


.0611 
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0.9 .6600 .4800 


.378 
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.0518 
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0.9 .8400 .3000 


.378 


.1805 


.1634 


.0171 
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.0174 



with g(p) = (sin^ ^ + Qcos^^). Inserting this /3 into the 
Boltzmann equation (4.1), we find that 

A = 0, B = y{vl + vl-l)g'{0). (4.7) 

As g'(6) is non-zero for a flattened model, /3 is an approxi- 
mate integral for those orbits for which the condition 

vl + vl^l, (4.8) 

is fulfilled. 

The easiest way to investigate the accuracy of the ap- 
porximate integral is to numerically integrate representative 
orbits and find the fluctuation in /3. The results of such 
computations are listed in Table 2, while Fig. 3 shows cross- 
sections of the orbits in the meridional plane. All the orbits 
are normalised to i? = and are identified by their starting 
values Ro and zo, with Vr and ve zero there. In each case, 
the maximum and minimum values of /3 along the orbit are 
recorded, together with the variation in /3 and in . The 
top section of this Table gives the data for four thin tubes in 
a mildly non-spherical potential (g = 0.9). For all these thin 
tubes — even for the last one with the highest ^-excursion 
— the fluctuations in /3 are small. The middle section of the 
Table contains two thin tubes with high ^-excursions, but 
for strongly flattened potentials (g = ^\/2 is the limit at 
which the density on the 0-axis vanishes). For the flattest 
potentials, the /3 excursions are becoming uncomfortably 
large. Finally, in the bottom section of the Table, three fat 
tubes — one very fat — show similarly large excursions in 
/3. The last column in Table 2 gives the variation in for 
the same nine orbits. Comparison of this column with the 
one for 8I3 shows that the excursions of /3 are persistently 
smaller than those of , in fact by a factor exceeding ten for 
all thin tubes with g = 0.9. We conclude that /3 is indeed a 
partial integral, in the sense that it is quite accurately con- 
served on thin and near-thin tube orbits in Binney's model 
with g > 0.8. It is not so well conserved for fat tubes. 

Table 2 contains a sequence of four orbits (B, E, F, and 
G) which all have E = and Lz = 0.378. They there- 
fore can not be discriminated by the classical integrals, even 
though they have distinctly different shapes, as shown in 
Figure 3, running from exactly thin (B) to very fat (G). But 
the approximate integral /3 can discriminate between them. 
Figure 4 shows the range in /3 (/3,min to /3,max of Table 2) 




R 



Figure 3. Cross sections of the orbits of Table 2 with a merid- 
ional plane {R,z) for the potential of Binney's oblate scale-free 
model with q = 0.9 at energy E = 0. The zero- velocity curve (for 
Lz = 0) is indicated by ZVC. Orbits A — D are thin tubes. Orbits 
E — G are fat (or thick) tubes that fill the indicated areas. The 
dot indicates the upper turning point of orbit K, which fills most 
of the area inside the ZVC (see section 5.2). 

for each of the four orbits B, E, F, and G. These ranges do 
not overlap. This shows that the average value of /3 for each 
orbit is a useful discriminant between orbits of the same E 
and Lz- The bottom half of Figure 4 indicates that L^ is a 
less effective discriminant than /3 , particularly between thin 
and near-thin orbits. 

4.3 Three-integral distribution functions 

The two exact integrals Li , L2 and the partial integral /3 
can be used to construct approximate DFs /(/i,/2,/3). We 
define a basic solution of the Boltzmann equation to be a 
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the density of our model can now be written as 
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Figure 4. The variations of and Is for orbits B, E, F, and G, 
which all have the same values of E and Lz. The average value of 
the approximate integral Is distinguishes all orbits quite clearly. 

is a much less useful discriminant. The exact range of Is and 

for each orbit is given in table 2. 

product of powers of the three integrals (cf. White 1985; 
Section 3.1) 



r j-k jn jm 

]k,n,m — Jl J2 



with 



1 



sin^" e 



j.2{k — n — m) g^~'^^0^ 



Fk 



(4.9) 



Fu,n,m = e-'''^""^""'^"lKf(vl + vlr. (4.10) 

The general solution of the Boltzmann equation can be writ- 
ten as a sum over the basic solutions with arbitrary coeffi- 
cients. We require the subset of the general solution consis- 
tent with the density of Binney's model, which is given by 
expression (2.6), with h(d) given in (2.9). To reproduce the 
r-dependence of this density, we must restrict ourselves to 
those basic solutions (4.9) with k = l-\-n-\-m. Furthermore, 
to reproduce the ^-dependence of the density, we must re- 
strict ourselves to just two values of n, namely ra = (i.e., 
k — m = l), and ra = 1 (i.e., k — m = 2). Implementing these 
restrictions and integrating over the velocities, we obtain 

1 v-^ r 1 c 4 ^^^^ ^ c 1 fA-i-i\ 

P=~il^ [^O.™;^'^".".™ + -^l.^-^^'^O.l.m I . (4.11) 



Here the ylo.m and A\^rn are arbitrary amplitudes and the 
So.o.m and 5o,i,m are the velocity integrals 



So n ; 



-Fl + n+m,n,md?)rd?)ed?),^, 



(4.12) 



which can be worked out explicitly (see Appendix B). Fi- 
nally, to make the density (4.11) implied by the phase-space 
DF exactly consistent with the density of Binney's model, 
we have to fulfill two conditions 



E 

m 

E 



2 - Q, 
2Q - 2. 



(4.13) 



If we replace the two sets of coefficients ylo.m and A\^rn by 
the weights FFo m and W\ m according to 



Ao,mSo,0,m — (2 — Q)Wo,m, 
Al,mSo,l,m = 2{Q - l)Wl,m, 

then the two conditions become 



Wo, 



1, 



(4.14) 



(4.15) 



f = J2{Ao,^ir^'lS 



I A Tm + 2 J jn 



(4.16) 



or, more explicitly, as 



g{0) So 



+ 2 



//O ia'^I.™ -(2+m)(f r+f^+fj) 2 / 2 , 2^r; 



(4.17) 



This DF should be used with caution. It involves Is which 
is a good partial integral for thin and near-thin tubes, but 
not for fat tubes. The solution (4.16) is exact in the limit 
in which it contains only terms with m = {Wo,m = and 
W^i.m = for m > 0). In this limit it depends only on the 
globally defined exact integrals Ii and I2 given in (4.5) and 
so reduces to the two-integral DF found earlier by Toomre 
(1982) and Evans (1993) through another approach. 

In the opposite limit of very large m, the DF (4.16) 
also seems trustworthy, for the following reason. For large 
m, the factor exp[— m(?)g -|- ?'|)](?'e + "I)™) which occurs in 
both terms of the DF (4.16), has its maximum at Vg -\- = 
1, and drops off sharply at both sides of this maximum. 
According to equation (4.8) this defines the part of phase 
space occupied by thin tubes. Thus, for large m, the solution 
(4.16) concentrates on the phase-space portion in which Is 
is an acceptably accurate integral. We conclude that the DF 
(4.16) may be used safely for any combination of m = and 
large m. 



4.4 Velocity dispersions 

The velocity dispersions which correspond to our solution 
for the phase-space DF can be found as follows. In (4.16), 
the velocity components occur only in their squares. Ac- 
cordingly, the off-diagonal terms of the velocity dispersion 
tensor vanish, i.e.. 



iVrVg) 



0, 



0, 



0. 



(4.18) 



The diagonal terms are obtained, as usual, by multiplying 
/ with v^j , (where j stands for r,d,4>, respectively) and in- 
tegrating over all velocities. We then find 



1 [ 2-Q ■S^J.O.m 

g(d) ^ 50,0, m 

m 

+ 2^(Q-i)ym,. ^"^■^ 



(4.19) 



So- 



where we have given the velocity integral (4.12) the more 
general form 



5, 



Vjll + r. 



jd?)rd?)ed?),^. 



(4.20) 



The velocity integrals can be computed analytically by in- 
troducing spherical coordinates in the velocity space. This 
transforms the three-dimensional integral into a product of 
three one-dimensional integrals (see Appendix B). The re- 
sults of these computations give for the three velocity dis- 



1_^2-Q 



Wo, 



g(e) 2(l + m 

m 



2\ 1 [2-Q , sin" e ,^ 



(4.21) 



Notice that and (?)|) are independent of the weights 

Wo,m and FFi,m, and hence are independent of m. At the 
poles, (vg) and (?)|) are equal, while at the equator, (v^) is 
larger than (vg). If the weights for m = are dominant (two- 
integral DF), the dispersions are isotropic at the poles, while 
at the equator (v^) and (vg) are equal, but (v^) is larger. 
If the weights for large m are dominant (predominance of 
thin tubes), then (v^) is much smaller than the other two. 
The weights FFo.m and FFi,m may be chosen independently 
of each other so that relatively small radial dispersions can 
be achieved on the poles or on the equator separately, within 
limits depending on Q. These models are similar to the Type 
1 tangentially anisotropic model 2 of Dehnen & Gerhard 
(1993). In the strongly flattened limit (Q = 2, g{0) = 1 + 
cos" 6), the first term in (4.21) vanishes for each dispersion. 
In the spherical limit (Q = l,g{9) = 1), the second term in 
(4.21) vanishes for each dispersion and we recover our earlier 
results (3.6) for the self-consistent isothermal sphere. 

5 COMPARISON OF JEANS AND BOLTZMANN 
APPROACHES 

5.1 Spherically aligned velocity ellipsoid 

The comparison of the two-parameter set of solutions (2.28) 
from the Jeans approach with the results (4.21) from the 
Boltzmann approach, shows immediately that the terms 
produced by the second parameter, H2, of the Jeans so- 
lution are not reproduced by our approximate Boltzmann 
solution. This is a direct consequence of the fact that our 
approximate partial integral contains only the squares of the 
individual velocity components, but not the cross term VrVe, 
just as is the case for the two classical integrals. Thus, any 
distribution function composed of these three integrals must 
produce a vanishing {vrVg), and spherically aligned velocity 
ellipsoids. We have already anticipated this circumstance by 
explicitly writing out the Jeans solutions for the case H2 = 
in (2.31). 

A comparison of equations (4.21) and (2.31) shows that 
the Jeans approach with H2 = — or F2{0) = — and the 
Boltzmann approach give identical results for (vg) and (v^), 
a reassuring circumstance. Now as to (v^), solutions (4.21) 
and (2.31) will give identical results for all 9, if the coeffi- 
cients of the g(d) terms and, separately, those of the sin" 8 
terms are identical. These two conditions can be brought 
into the form 



1 + m 

m 

^ 2 + m 



1 + 



Q{2-Q}' 
Hi 



(5.1) 



mere are many ways 10 saiisiy inese consirainis. ror ex- 
ample, let us choose in each sum only two non-vanishing 
terms, namely one with m = (two-integral solution) and 
the other with m ^ 00 (thin tube solution). Then 



Win 



«0, 
"1, 



Wo, 
Wi, 



■ «o, 

■ «i. 



phich gives for the two sums 



Wo 



= «0, 



2 



Wi 



^ 2-F 



■ = a\. 



(5.2) 



(5.3) 



Conditions (5.1) then determine ao and a\ for any given 
R\ . However ao and a\ will fall into their physical range (0 
to 1) only if 



-\Q{2-Q) <Hi <0. 



(5.4) 



The lower limit for Hi given by condition (5.4) is fulfilled by 
the Jeans solution according to equation (2.34). However, 
the upper limit for Hi given by (5.4) is not required by the 
Jeans solution and arises from our Boltzmann approach, in 
which we took care of the thin tubes, but not of the fatter 
ones. 

We conclude that the Boltzmann approach with our 
partial integral reproduces the two-parameter solution of 
the Jeans equations accurately, but only if H2 = 0, and 
if Hi < 0. It would seem plausible that these limitations 
could be raised by, a), developing our third integral to con- 
tain a VrVe term, and by, b), finding another partial third 
integral that should be of good accuracy for flat tube orbits, 
with moderate ^-amplitudes. These two extensions are not 
attempted here. 



5.2 The cross-term 



The general solution (2.17) of the Jeans equations for oblate 
scale-free potentials contains many solutions with F2{0) ^ 0, 
i.e., (vrVe) ^ 0, which have velocity ellipsoids that are not 
aligned with the spherical coordinates r and 9. Such a term 
arises when the DF depends on a third integral which con- 
tains a term proportional to VrVe. Might a more accurate 
partial third integral than our I3 = g{d) contain such a 
term? We have not found the overall answer to this question. 
The early numerical experiments by Richstone (1984) for the 
oblate scale-free logarithmic potential suggest that many 
self-consistent solutions show misalignment of the velocity 
ellipsoid relative to {r,d). Levison & Richstone (1985a, b) 
found the same for r~'^ densities in Binney's potential. Here 
we add a simple argument indicating how such a misalign- 
ment might come about. 

Consider first the solution constructed exclusively from 
thin tubes. Such solutions have been derived in detail for 
Stackel potentials (Bishop 1987; de Zeeuw & Hunter 1990). 
A thin tube has a meridional cross section which is elongated 
— larger r — in the polar direction of the oblate potential 
(see orbits A, B, C, and D in Fig. 3). The velocity vector 
{vr, Vg) must fall along this cross section. The polar elonga- 
tion then gives Vr-vg < for all points above the equator (the 
opposite sign holds for points below the equator). Thus for 
the thin-tube solution (vrVg) < above the equator. Note 
that this sign for (vrVg) holds for the scale-free logarithmic 
potential and is therefore not caused by the existence of a 
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Figure 5. Contributions to (vrVg) by individual orbits as a func- 
tion of the angle 6. The near-thin orbits A, B, C, D all contribute 
negative values for all 6. They are shown in a broken line. Very 
fat orbits, such as F,G and K, contribute both positive and neg- 
ative values. They are shown in an unbroken line. 

regular core, such as occurs in Stackel models. 

In contrast, the two-integral solution, which contains 
all orbits, not just the thin tubes, has {vrVg) = 0. Do the 
thicker tubes then contribute positively to {vrVg) above the 
equator? The answer seems to be yes, as indicated by the 
following evidence. 

A solution for a DF in a given potential can be thought 
of as consisting of the sum of contributions from all the in- 
dividual orbits involved. In an oblate potential a tube orbit 
provides a {vrVg) contribution that is a function of r and 9. 
In a scale-free oblate potential, for a scale-free DF, a spe- 
cific orbit can be taken to represent a sequence of similar 
orbits, all equal except in scale, i.e., in energy (Richstone 
1980; Schwarzschild 1993). This whole sequence is taken 
into account by projecting the integrated orbit onto a ref- 
erence circle in the meridional plane, with the appropriate 
scaling rules. For the logarithmic scale-free potential the 
rules are: Leave the velocities unchanged and weigh each 
point in the integrated orbit by dt/r. Thus the resulting 
(vrVe) for each orbit, i.e., orbit sequence, depends only on 
0. 

Results of (vrVe) computations are given in Figure 5 
for the seven orbits of Table 2 with q = 0.9. The four 
thin tubes A, B, C, and D give negative contributions to 
(vrVe) for the entire range of 9, just as expected from the 
geometrical argument given above. The absolute values of 
(vrVe) for these thin tubes are relatively small, generally 
less than 10% of (vg), leading to less than 6° of tilt for the 
velocity ellipsoid in thin-tube solutions. For lower q values 
(galaxies flatter than E3) the maxima of K^r^e)! are higher, 
reaching about 20% for q = 0.8, and still larger values for 
smaller q, such as used by Richstone (1984) and Levison & 
Richstone (1985a). 
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thinnest of the three, is still thin enough to give all nega- 
tive (vrVe), like the thin tubes. Orbit F shows in Figure 5 
the first positive {vrVg), but still very small. Only orbit G 
of all shown in Figure 3 makes substantial positive contri- 
butions to (vrVe), and that mainly near its upper margin 
(near the smallest 9 it reaches). To show that orbit G is not 
a rare case, we have added a very thick orbit K (its upper 
turning point is plotted in Figure 3). Figure 5 shows that 
orbit K makes strong positive contributions to {vrVg), again 
concentrated near its upper margin. This shows explicitly 
that thick orbits like G and K balance the thin orbits to give 
(vrVe) = for the two-integral solutions. 

Can we understand the positive contributions of the 
thick orbits in simple geometrical terms? The direction of 
the upper margins of tube orbits in the meridional cross sec- 
tion shown in Figure 3 gives the answer. These upper mar- 
gins are not exaclty radially directed but are slightly turned 
clockwise — reminiscent of the ellipsoidal coordinates that 
define the margins of orbits in a Stackel potential. This 
tilt in direction gives a positive {vrVg) on the upper margin. 
In the neighbourhood of this margin the velocity vector of a 
thick tube orbit straddles the direction of the margin. Hence 
near the upper margin a thick tube orbit should contribute 
a positive {vrVg), just as Figure 5 shows for orbits G and K. 

We conclude that for the oblate scale-free logarithmic 
potential: A solution made up of thin or near-thin tubes — 
tangentially anisotropic — will have a negative {vrVg) above 
the equatorial plane, though of modest absolute amount. 
In contrast, a solution made up of thick tubes — radially 
anisotropic — may have positive {vrVg) above the equator, 
at least in part of the meridional plane. 



6 OBSERVABLES 

6.1 General scale-free case 

We define a coordinate system (x' , y' , z') with z' along the 
line of sight by the transformation 

/x'\ / 1 \ / x\ 

I J/' I = I — cos i sin i I I ?/ I • i^-^) 
\ z' J \ sin i cos i J \ z J 

Here i is the inclination, defined such that i = 7r/2 corre- 
sponds to edge-on observation. Integration of model proper- 
ties over z' then gives the projected (observable) properties 
on the sky-plane (x',y'). General expressions for the line- 
of-sight projected properties can be found in Appendix A of 
Evans & de Zeeuw (1994). 

The projected properties of scale-free models with po- 
tentials (2.5) and densities (2.6) are themselves scale-free. 
It is therefore useful to define polar coordinates (R',Q') in 
the plane of the sky by 

x' = -R'sm&', y' = R'cos&', (6.2) 

so that R'^ = x'^ + y'^ and 9' runs counterclockwise, and is 
measured from the j/'-axis ("North"). Then the projected 
surface density is 

S(j?',e')= (6.3) 
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second moments. 
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o„ (s sin i — COS 0' COS i) +sin 0' 
tan 61= ^ ' 



ds h(0[Q',s]) 
S-' gH0[e',s])' 



(s cos i + cos 0' sin i)2 



(6.4) 



(6.5) 



Similar expressions for the projected velocity moments can 
be written down using equations in Appendix A of Evans & 
de Zeeuw (1994). 



6.2 Binney's model 

The motivation for introducing the two-parameter Jeans so- 
lution (2.28) for Binney's model in Section 2.6 is that the 
line-of-sight integrations can be done exactly to give the ob- 
servables. We have shown that a one-parameter subset of 
these Jeans solutions {H2 = and Hi < 0) corresponds to 
physical DFs, at least to good approximation. Let us now 
see how these galaxy models appear to an observer. 

The projected surface density of Binney's model is 



E(x', y') 



g'3 [x'2 -I- j/'2g/-2]3/2 

7rg 1 



(6.6) 



g'3i?' (cos2 0' + g'-2 sin2 0')3/2 ' 



with g'2 = cos2 i -\- srn^ i. The six components of the 
projected second moments follow from equation (A7) of 
Evans & de Zeeuw (1994). The line-of-sight second moment 
(v^i) = (j'fos) = (n^^i^i) is the best known of these quantities. 
The other two diagonal components and (n^iyi) are 

the dispersions in the proper motions. The off-diagonal com- 
ponents, (fi%yi), and ^i) measure the correlations 
between the distributions of proper and radial velocities. All 
these quantities are of the form 

2 ^ aijx'2 -I- b,jy''^ + Ctj x'y' 



4(x'2-Kj/'2) ' (6.7) 

cos^ 0' -\-hij sin^ 0' -\-Cij cos 0' sin 0'), 



where the non- vanishing ajj, htj, and Cjj depend on i, q and 
Hi and are given in Table 3. 

The velocity profile (VP) is the distribution of line-of- 
sight velocities vios = v^i measured at position (x' , y' ) on the 
sky (see, e.g., van der Marel & Franx 1993; Gerhard 1993). 
It is obtained by integrating the DE along the line-of-sight 
and over the tangential velocity components (see e.g., Evans 
& de Zeeuw 1994). The normalised VP is 



VP(v^i,x' , y')-- 



1 



E(x',j/') 



dV,, dVy,f(h,h,h).((<.^) 



Recently, there have been a number of investigations of the 
range of shapes of the VPs of flattened galaxy models with 
two-integral DEs (e.g., Evans 1994; Qian et al. 1995; CaroUo 
et al. 1995). The variety of shapes of the VPs of three- 
integral DEs is expected to be richer still (e.g., Dehnen & 
Gerhard 1993). 

Eirst, we consider the VPs of the DE (4.16). When only 
the m = terms are present, it reduces to the two-integral 
DE. In the non-rotating case, the VP of the two-integral 
DE is the sum of two Gaussians, as shown in Appendix B 
of Evans & de Zeeuw (1994). It is useful to contrast this 
with the VP of a model dominated by near-thin tube orbits, 
say with just the m = 20 terms in (4.16). Eig. 6a shows 
the VPs observed on the major axis for edge-on viewing 
(9 = 90°) in these two cases. The intrinsic flattening of 
the model is E3 (g = 0.9). When this galaxy is built from 
near-thin tubes, the VP becomes squatter. The wings of 
the VP become shallower as there are fewer eccentric orbits, 
while there is a sudden and substantial rise in the VP at 
the circular speed. In the commonly-used Gauss-Hermite 
decomposition (see e.g., van der Marel & Eranx 1993), the 
VPs of the near-thin tube models possess negative values 
for the hi coefficient. The dotted curve shows the VP of the 
m = 1 component DE, which has properties intermediate 
between the two-integral and near-thin tube DEs. 

Second, we investigate the VPs of the DEs correspond- 
ing to our Jeans solutions. We choose the simplest con- 
ditions (5.2). Eor any flattening, there is a one-parameter 
family of Jeans solutions depending on Hi . Eor any Jeans 
solution, there is a one-parameter sequence of DEs, depend- 




Figure 6. (a) VPs of Binney's model with an ellipticity of E3 viewed edge-on [i = 90°). The viewing point is on the projected major 
axis. The full curve is the VP generated by the two-integral DP (m = in (4.16)). The dashed curve is the VP of a model dominated 
by near-thin tubes (m = 20). The dotted curve shows an intermediate case (m = 1). (b) VPs of Binney's model which are consistent 
with the Jeans solutions (2.31) with Hi = —0.45. The full curve is built from the m = and m = 1 components, the dashed curve from 
the m = and m = 20 components. 



ing on m. In other word, the DF includes a contribution 
from the two-integral DF (m = component), together with 
a contribution from a component at one other m. Fig. 6b 
shows the VPs generated when the second component is 
m = 1 (full curve) and m = 20 (dashed curve). Both 
VPs have the same second moments and correspond to the 
same Jeans solution (2.31), which is tangentially anisotropic 
(Hi = —0.45). The effect of including lots of thin and near- 
thin tubes in the DFs is to produce characteristic kinks in 
the VP at roughly the circular speed. An observer sees sub- 
stantial numbers of thin and near-thin tubes with v,f, Si ±1 
and with ve Si ±1. When the models are viewed edge-on, 
this produces kinks in the VP at v^i Si ±1 as well as the 
peak at v^i Si 0. When the models are viewed pole-on, the 
observer sees mainly the ve velocity components of orbits 
at the tangent point. So, the characteristic appearance of 
the VP of a model dominated by near-thin tubes is robust 
under inclination angle changes. 

All the VPs of the models with many near-thin tubes 
have negative hi coefficients in the Gauss-Hermite decom- 
position. This result is valid for a wider range of flattened 
scale-free models. For example, CaroUo et al. (1995) showed 
that the VPs of f{E, Lz) non-rotating scale-free spheroids 
with a range of density profile slopes have negative hi coef- 
ficients. Increasing the tangential anisotropy by populating 
more near-thin tube orbits at the expense of the fat tubes 
will lead to /t4-values that are more negative. 



7 CONCLUSIONS 

We have derived the general solution of the Jeans equations 
for arbitrary oblate scale-free densities in oblate scale-free 



logarithmic potentials. This provides all possible second mo- 
ments (vr), (vrVe), (vg) and (v^) that can support the given 
density in the gravitational field generated by the potential. 
The general solution contains two free functions, Fi{9) and 
F2(0), which give the angular variation of (?)^) and {vrVe), 
respectively. We have studied the specific case of Binney's 
model in detail. This is the self-consistent scale-free model 
with spheroidal equipotentials that is often used as a simple 
dark halo model. We have identified a simple two-parameter 
set of Jeans solutions for this model with the useful prop- 
erty that the observable kinematic quantities can be given 
explicitly. 

In order to identify which of these Jeans solutions cor- 
respond to physical phase-space distribution functions, we 
first investigated the spherical limit. In this case the poten- 
tial reduces to that of the singular isothermal sphere, and 
all physical Jeans solutions for scale-free densities of arbi- 
trary flattening must have F2{0) = 0, i.e., {vrVg) = 0, so 
that the velocity ellipsoids must be spherically aligned. We 
have provided DFs for all the remaining Jeans solutions. 

We then investigated the problem of finding distribu- 
tion functions associated with the Jeans solutions in the 
flattened scale-free logarithmic potentials of Binney's model. 
We found an approximate solution of the collisionless Boltz- 
mann equation, and showed by numerical orbit integration 
that it provides a third (partial) integral of good accuracy 
for thin and near-thin tube orbits. It is a modification of 
the total angular momentum. We then constructed sim- 
ple three-integral distribution functions, and showed that 
the kinematic properties of these approximate DFs agree 
with the subset of our two-parameter Jeans solutions which 
have spherically aligned velocity ellipsoids (i.e., the param- 
eter H2{0) = 0). Although our DFs are approximate, this 
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tened anisotropic scale-free Binney models can be approxi- 
mated quite accurately by assuming spherical alignment of 
the velocity ellipsoid. 

The derivation presented in Section 2 suggests that the 
general solution of the higher-order moment equations may 
be found in a similar way. We expect that at each order fur- 
ther free functions of 9 appear. The problem of picking those 
that correspond to physical DFs remains. One might also 
attempt to solve the entire coUisionless Boltzmann equation 
directly. White (1985) was able to find the general solu- 
tion for scale-free spherical models because in this geometry 
three exact integrals of motion are known. Since there is no 
exact third integral for an arbitrary axisymmetric scale-free 
potential, an exact solution of the coUisionless Boltzmann 
equation seems to be precluded, even though it may be pos- 
sible to find solutions of the moment equations to arbitrary 
high order. It is still possible that exact, globally defined 
third integrals may exist for some specific choices of g(d). 

Many of the properties of the two-integral Binney model 
carry over to the general family of axisymmetric power-law 
galaxies (Evans 1994), which have potentials of the form 



$ : 



l3{Rl + R^ + ; 



11 ) 



for /? = 0, 



for /? / 0. 



(7.1) 



Compared with the case considered in this paper, these mod- 
els have two extra parameters: the power /? of the radial 
dependence of the spheroidal potential, and the core radius 
Rc- The general solution of the Jeans equations given here 
for Binney's model (Rc = 0,/? = 0) can be generalized to 
scale-free cases with /? 7^ 0. When Rc ^ d the radial de- 
pendence of the stresses is not known a priori, and solving 
the Jeans equations becomes a harder problem. Our results 
suggest that assuming spherical alignment of the velocity 
ellipsoid might be a good approximation over most of the 
model, and the general solution with {vrVg) = is available 
(Bacon 1985). The two-parameter set of solutions presented 
in Section 2.6 can be generalized to the /? 7^ 0, i?c 7^ case, 
and the observable properties of these solutions can again 
be written down explicitly. The partial integral presented 
here for Binney's model can also be modified to the gen- 
eral power-law models. These results will be discussed in a 
future paper. 

A further generalization is to consider triaxial models. 
The axisymmetric power-law galaxies can be made triaxial 
by taking the equipotential surfaces to be similar triaxial 
ellipsoids rather than spheroids. There are then three Jeans 
equations for the six components of the stress tensor. The 
scale-free case can be treated analogously as we have done 
here, and the two-parameter solution can be generalized as 
well. We expect it to contain three free parameters. Analysis 
of this interesting case is in progress. 
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APPENDIX A: GENERAL SOLUTION IN 
CYLINDRICAL COORDINATES 

In many studies, properties of axisymmetric dynamical mod- 
els are given in terms of standard cylindrical coordinates 
{R, 4>, z). Our general solution (2.17) results in the following 
components of the stress tensor in cylindrical coordinates: 



p{vr} 



T-'gHe)' 



p{vRVz) 



G2{0) 
r'lg^iey 

T-ig^{ey 



(Al) 



Willi / = TL -t" -S , ^-^4 = -^4 aiiu 



Gi{0) = _Fi(6i)sin^6i + 2_F2(6')sin6icos6i + _F3(6')cos^6', 
G2 (61) = [Fi {0)-F3 (61)] sin cos 6 

+ _F2(6)[cos^6-sin^6], 
63(6) = _Fi(6)cos^6-2_F2(6)sin6cos6 + _F3(6)sin^6. 

Spherically aligned solutions have F2{0) = 0. Cylindrically 
aligned solutions, (vrvz) = 0, must have 62(6) = 0, and 
hence are obtained by taking Fi{0) = I{0) + _F2(6)tan6, 
with F2{0) arbitrary. Then 

«^W=^W + ;i5£6' «3(6) = /(6). (A3) 

Taking also F2{0) = gives the / = f{E, Lz) solution with 
(vji) = (vl) and {vrv^} = 0. 



APPENDIX B: AN AUXILIARY INTEGRAL 

The normalising constant in our three-integral DF involves 
the integrals (4.12) : 

So,n,m = dvrdvedvd, 

J J J (Bl) 

^_(l + „+„)(,2 + ,2 + ,J)^2n^^2 ^ ^2^™^ 

To evaluate, we introduce spherical polars in velocity space 
(see BT, p. 293) 

vr = vcosr], ve = v sin rj costp, v ,f, = v sin rj sin tp . (B2) 

This reduces the three-dimensional integral to a product of 
three one-dimensional integrals, which can be evaluated us- 
ing formulae 3.621.5 and 3.381.4 of Gradshteyn & Ryzhik 
(1980). We finally obtain: 

7rr(ra-|-i)r(ra-|-m-|-l) , , 

JO.n.m — 5 . i-tS-jj 

The generalised velocity integrals (4.20) can be found in a 
similar manner. We find 
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